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ARE GALAXY PROPERTIES ONLY DETERMINED BY THE DARK MATTER HALO MASS? 
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ABSTRACT 

Conditional luminosity function (CLF) of dark matter halos supersedes simple models on the number 
of galaxies as a function of the halo mass, or the so-called halo occupation number, by assigning a 
luminosity distribution to those galaxies. Suggestions have now been made that both the clustering 
strength and properties of hosted galaxies of a given dark matter halo depends also on the age and the 
environment of that halo in addition to halo mass. Based on simple CLF models, we find that at least 
one suggestion that made use of central galaxy properties may be affected by uncertainties in the group 
catalog from which clustering properties were measured to address the age dependence. Establishing 
how galaxy properties of fixed mass halos change with the environment or the age is challenging, if not 
impossible, given uncertainties in group and cluster catalogs. It may be possible to address this issue 
through statistics based on a combined study of luminosity and color distributions, and luminosity- and 
color-dependent galaxy clustering, all of them as a function of the galaxy overdensity, though no single 
statistic is likely to provide the ultimate answer. The suggestion that the age dependence of the halo 
bias invalidates the analytical models of the galaxy distribution is premature given the improvements 
associated with modeling approaches based on conditional functions. 

Subject headings: 



1. INTRODUCTION 

Recently, several claims have been made that the galaxy 
properties of a given dark matter halo may depend on the 
environment in which the halo resides (Gao et al. 2005; 
Harker et al. 2005; Yang et al. 2005b). A connection is 
generally made to the age of the halo since age may be the 
primary factor that determines the environment in which 
the halo forms. Numerical simulations of the dark matter 
distribution suggests that older halos are more strongly 
clustered than younger halos of a given mass (Gao et al. 
2005). Such an age dependence on the halo bias is not yet 
present within successful analytical descriptions of the halo 
distribution and clustering properties (e.g., Mo & White 
1996; Sheth, Mo & Tormen 2001). These results invalidate 
analytical models of the galaxy distribution since these 
models assume that the galaxy properties are only deter- 
mined by the halo mass (Kravtsov et al. 2004; Cooray & 
Sheth 2002 for a review). 

Beyond simulations, whether galaxy properties of a 
given halo mass vary with the environment is yet to be es- 
tablished observationally. The galaxy morphology-density 
(Dressier 1980; Dressier et al. 1997) and the color-density 
(Blanton et al. 2004) relations do indicate that galaxy 
properties change with the environment: red, early-type, 
galaxies are usually found in denser environments than the 
environments where blue, late-type, galaxies are found. It 
could, however, be that the increasing fraction of early- 
type galaxies with increasing density is simply a reflection 
of a higher fraction of early-type galaxies in more massive 
halos (e.g., Cooray 2005a; Weinmann et al. 2005). Thus, 
interpretations of the density-morphology or density-color 
relations are generally inconclusive with some suggesting 
support for the halo model (Blanton et al. 2004) while 
others not (Goto et al. 2003; Balogh et al. 2004). 

To study the difference between galaxy properties of 
same mass halos when either the environment or the age is 



varied one must find fixed mass halos while a second vari- 
able, such as the environment, is varied. The test related 
to the environment could be achieved with large catalogs 
of galaxy groups, but since one may not be able to reli- 
ably assign a mass to an individual group, it is unlikely 
that there is an easy way to uniquely identify whether one 
is comparing halos of similar mass. Since the expected dif- 
ferences are likely to be at the few to, at most, ten percent 
level (see e.g., discussion in Zheng & Weinberg 2005), the 
mass assignments must be accurate to the same precision. 
Statistics that measure properties after randomizing galax- 
ies can partly address the question on the environmental 
dependence by testing whether the surrounding large-scale 
environment is responsible for properties within a halo. It 
could be that the internal structure of a halo, such as the 
dark matter concentration parameter, changes with the 
environment and galaxy properties reflect those changes. 

The second approach to address this question involves 
statistical study of a large sample of galaxies, such as clus- 
tering, binned in terms of an estimator of age for galax- 
ies. Using a catalog of galaxy groups and clusters from 
the 2-degree Field Galaxy Redshift Survey (2dFGRS; Col- 
less et al. 2001), Yang et al. (2005b) recently considered 
the cross-correlation function between the central bright 
galaxies in galaxy groups and the general population of 
galaxies above a certain luminosity. For groups with halo 
masses above Af*, they found the groups with a lower lumi- 
nosity central galaxy to be more strongly clustered than 
the ones with a higher luminosity galaxy. Furthermore, 
galaxies that are more passive in a given mass bin were 
found to be more correlated. Thus, a claim was made that 
the cross-clustering measurement in Yang et al. (2005b) 
indicates age and environmental dependence of halo bias. 

Using models of the conditional luminosity function 
(CLF; Yang et al. 2003, 2005a; Cooray & Milosavljevic 
2005b), we make a simple model of the cross-correlation 
between central galaxies of groups and the galaxy distri- 
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bution and suggest that the measurement in Yang et al. 
(2005b) may be contaminated by uncertainties in their 
group catalog involving mass assignments. We find that, 
for example, if one selects roughly 10% to 20% fraction 
of satellite galaxies of larger mass halos mistakenly as the 
central galaxy in lower mass halos, then the bias factor for 
smaller halos with fainter luminosities is boosted. With 
a simple relation between the starformation rate and the 
galaxy luminosity, the clustering results binned in terms 
of starformation rate may also be explained without re- 
quiring a second parameter. It is unlikely that Yang et 
al. (2005b) measurements provide a definitive answer on 
the age dependence of galaxy properties of fixed halo mass 
with age captured by the either the star-formation rate or 
the luminosity of the central galaxy. 

Here, we also consider if there is an appropriate sta- 
tistical test on the density and age dependence of galaxy 
properties within a given halo that is affected less by sys- 
tematics. While we have yet to come up with a simple test 
to provide the ultimate answer a combination of measure- 
ments on the luminosity and color distributions, combined 
with clustering, all as a function of galaxy density may be 
utilized to address this issue. Any single statistical mea- 
surement is affected by two potential variations: the mass 
function that varies the halo distribution with the envi- 
ronment and galaxy properties within a fixed halo mass 
that change with environment: 
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where both the mass function, dn/dM, and the CLF, 
$(_L|M), varies as a function of the environment cap- 
tured by (5 ga i, which is defined in terms of the density of 
galaxies in some fixed volume. To address the question 
whether $(L|M, <5 ga i) = $(L|M) observationally, the sta- 
tistical methods must allow a separation of effects coming 
from dn/dM(8 ga \) with <J>(L|M, <5 ga i) (see, e.g., Mo et al. 
2004; Cooray 2005a). This degeneracy is also what affects 
interpretation of the density-morphology or density-color 
relations. 

The Letter is organized as following: In the next section, 
we briefly summarize a model for the cross-clustering be- 
tween central galaxies in galaxy groups and the general 
galaxy distribution. We compare our model predictions 
with those of Yang et al. (2005b) and suggest that their 
observations can be explained in terms of a small contami- 
nation in the group catalog. Furthermore, we comment on 
potential tests to address whether <!>(L|M, <5 ga i) = <1>(L|M) 
or not. Finally, as a side note, even if it turns out that 
the CLF has a second parameter, we will argue that this 
does not imply that the analytical halo model is wrong, 
but rather it can be easily improved by including the envi- 
ronmental dependence when describing galaxy properties 
through conditional functions. 

2. GALAXY BIAS OF CENTRAL GALAXIES IN GROUPS 

In Yang et al. (2005b), the cross-correlation between 
luminous galaxies (presumed to be at halo centers) and a 
large sample of galaxies, above some luminosity, is studied 
as a function of the luminosity and star-formation rate of 
the central galaxy. In models based on the conditional 
luminosity function, the bias factor of galaxies above a 



minimum luminosity L m ; n can be calculated as 
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where the galaxy LF is $(£) = j dM dn/dM $(Z,|M) 
and bha,i (M) is the halo bias. We make use of Sheth & 
Tormen (1999) description for the mass function and the 
halo bias, consistent with numerical simulations for both 
the halo mass (Jenkins et al. 2001) and average bias (Sel- 
jak & Warren 2004). This CLF $(L\M) can be divided 
to central and satellite galaxies following the approach ad- 
vocated by Cooray & Milosavljevic (2005b). The central 
galaxy luminosities are assigned to halo masses through 
the Cooray & Milosavljevic (2005a) relation as appropri- 
ate for the 2dF bj band (Cooray 2005a; also, Vale & Os- 
trikcr 2004) including a scatter with a log-normal disper- 
sion of 0.17 as measured from SDSS LF and clustering 
data (Cooray 2005b). 

We write the central galaxy bias when divided to mass 
bins as a function of the central galaxy luminosity or star- 
formation rate rj as 



b c (L or r/| M min , M max ) = 
/ dM dn/dM &haio(M) [h{M, M, 
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where h(M, M m i n , M max ) is a function such that h(M, M m - m 
M max ) = 1 when M min £ M £ M max . While the 
group catalog is divided to bins in estimated halo masses, 
we allow for an uncertainty in both the minimum and 
maximum mass of each bin since it could be that halos 
with true masses below and above are assigned incorrect 
masses within the bin through statistical uncertainties. 
We achieve this with an analytical function that is essen- 
tially flat across the mass range with tails both above and 
below the mass cut-offs, and normalized to be unity. In 
addition to mass limits, we also allow for a slight varia- 
tion in the luminosity cut-off. In Yang et al. (2005b), the 
galaxy sample used involve a volume limited sample with 
galaxies M < — 19.45 + 5 log h in the &j-band. In our CLF 
models, instead of a sharp cut-off at L m i n corresponding to 
this absolute magnitude cut, we allow for a slight scatter 
in this luminosity as well. 

The predicted large-scale clustering is shown in Fig. 1. 
Here, we plot the product of b g (> L min )b c (L c \M min , M max ) 
as a function of L c , the luminosity of central galaxies, for 
mass bins defined by M m ; n < M < M max . Instead of all 
mass bins considered in Yang et al. (2005b), we consider 
two bins at the low and high mass end. The solid and long- 
dashed lines show the expected bias factor for two bins 
between 10 12 to 10 12 5 and 10 13 5 to 10 14 M . When com- 
pared to the measurements from Yang et al. (2005b), CLF 
models based on standard description for halo bias and 
mass function show a lower bias factor for low luminosity 
galaxies assigned to higher mass groups. Similarly, we also 
consider measurements divided based on star-formation 
rate (rj) . Since our model is based on luminosities we need 
a mapping between L c and n c . We considered a simple one 
to relation of the form n c = a(L c /10 10 Lq) + j3 and found 
out that combinations of (—5,1.5) and (—8,3.5) describe 
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Fig. 1. — The product of bias factor of central galaxies in groups b c (L c ) and galaxy bias for a volume limited sample above a bj magnitude 
of —19.45 + 51og?i with measurements subdivided to mass bins in galaxy groups that host the central galaxy. The measurements arc from 
Yang et al. (2005b) and we show their data for two mass bins from 10 12 to 10 125 (open symbols) and 10 135 to 10 14 Mq (filled symbols). In 
the left panel we consider measurements as a function of central galaxy luminosity while in the right panel we consider their measurements 
as a function of star-formation rate. The solid and long-dashed lines show the expected bias factor for central galaxies in these mass bins 
as a function of the central galaxy luminosity. The dot-dashed line shows a simple modification to the bias factor calculation where we 
include a small correction to central galaxy sample by allowing satellites in larger mass halos (> 10 14 Mq) to be identified (wrongly) as 
central galaxies and be assigned to a group with a mass in the range between 10 135 Mq to 10 14 Mq. The contaminant fraction / (see text) 
needed to explain left and right panels differ by a factor of 2 and could indicate extra scatter in mapping between luminosity, which is the 
primary ingredient in the CLF-based models, to the star-formation rate of central galaxies. The simple analytical description, which does not 
make any modifications to the halo bias as a function of age, suggests that these measurements with the observed galaxy distribution do not 
conclusively establish an age dependence to galaxy clustering properties. 



the data best in the high- and low-mass bins respectively. 
The difference in the mapping between L c and r\ c for the 
two mass bins suggest that starformation rate of a central 
galaxy depends on the mass of the halo in which it resides. 

To study the bias factor in the case of contaminations, 
we modify square brackets in Eq. 3 as [h(M, M m ; n , M max ) 
$ C (L|M)] -> [/i(M,M min ,M max ) * c (i|M) +/$ a (L|M) 
6(M, Mi, M 2 )], where 6(M, Mi, M 2 ) is the unit step func- 
tion such that either 9(M, Mi, M 2 ) = 1 when Mi > M > 
M 2 or zero otherwise. The assumption here is that some 
fraction of central galaxies are satellites of larger mass ha- 
los (above the mass range corresponding to a given bin), 
but mistakenly taken to be central galaxies by the group 
finding algorithm, say, and assigned with a lower halo 
mass. With Mi = 10 14 and M 2 = 10 14 5 , / - 0.002 
leads to the dot-dashed line shown in Figure 1 (right panel), 
while a value of / ~ 0.001 describes the dot-dashed line in 
Figure l(left panel) and with L c to r\ c mapping described 
above. The difference in / between the two plots may be 
suggestive of extra scatter in the mapping between L c and 
rj c . Note that / should not be interpreted simply as the 
number of central galaxies that may be contaminated in 
the group catalog since in our model description, a large 
difference exists between CLFs for central and satellite 
galaxies. Integrating over the LF, the difference in terms 
of a fraction of galaxies is at the level of 12% percent in 
the lowest luminosity end of the mass bin between 10 13 ' 5 
Mq to 10 14 Mq but drops to the level of a few percent 
and below as the luminosity is increased; this luminosity- 
dependent change in the contamination is not a reflection 
of a complex situation but rather the fact that in the mass 
bin 10 13 ' 5 Mq to 10 14 Mq there are more "true" central 



galaxies at the high luminosity end. 

It is unlikely that due to a few to ten percent uncertain- 
ties in any group catalog, based on statistical errors and 
systematics, one can use them directly to identify varia- 
tions in clustering properties between groups of different 
ages or environments. It is unclear if the measurements in 
Yang et al. (2005b) can be interpreted as a conclusive test 
on the age dependence of the galaxy or halo bias. One can 
question why the same situation, or increase in bias, does 
not happen in the low mass bin, for e.g., between 10 12 to 
10 12 ' 5 Mq as shown in Fig. 1. We find that if the contam- 
ination comes from adjacent mass bins, when M < M*, 
due to the fact that the halo bias is more or less constant 
at the low end (e.g., Mo & White 1996; Seljak & Warren 
2004), there is no significant difference in bias from the 
contaminants. Small changes are more visible in the high 
mass end both due to small density of galaxies and the 
large changes to the halo bias factor with a small change 
in the halo mass. It will be interesting to see if these vari- 
ations will also show up in similar tests with dark matter 
catalogs directly from numerical simulations whose mea- 
surements of dark matter properties will also be affected 
by small uncertainties. 

3. IS THERE A MORE RELIABLE TEST ON THE AGE 
DEPENDENCE OF GALAXY PROPERTIES? 

Instead of addressing whether bias itself is age depen- 
dent in the observed data, perhaps, it may be easier to 
address if galaxy properties of a given dark matter halo 
are defined only by the halo mass regardless of the age 
or the environment of the halo. Analytically, the test is 
simple: Is $(L\M,X) = $(L|M), or is there a second pa- 
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rameter X that determines galaxy properties of a halo in 
addition to the mass of that halo? This second parameter 
may be the environment, characterized by the overdensity 
Jgai, age at which the halo formed t/, or another param- 
eter that may determine galaxy formation and evolution 
(such as heating associated with reionization) . It is now 
well known that the satellite LF of groups and clusters 
are different at the faint-end; one finds a larger fraction of 
dwarf galaxies in clusters while the faint-end LF of groups 
is flat. Whether the LF is different in clusters and groups 
when the age or the environment of that group or cluster is 
varied may provide a useful test on the extent to which the 
assumption that the CLF depends only on mass is valid. 

Similarly, it may be useful to study differences in prop- 
erties of "fossil" systems relative to all other groups of the 
same mass. Fossil systems are expected to be older, iso- 
lated, and more biased relative to groups of the same mass. 
Perhaps the difference in properties between fossil and all 
other groups capture the extent to which $(L\M,tf) de- 
pends on the age of the halo t f . Since the fraction of fossil 
systems at group mass scales are ~ 10% (Milosavljevic et 
al. 2005), one should only expect variations in the galaxy 
sample roughly around the same level. 

One can easily address the issue on how large-scale en- 
vironment affect galaxy properties by randomizing the 
galaxy distribution within halos and comparing statis- 
tics before and after the redistribution. This was con- 
sidered by Yoo et al. (2005) using a simulation and a 
semi-analytical model and found <J 5% differences in clus- 
tering statistics involving the galaxy-galaxy and galaxy- 
mass cross-correlation functions. Whether these changes 
are a real reflection of galaxy properties and evidence for 
$(L|M, 5g a i) rather than $>(L\M) or a reflection of uncer- 
tainties in simulations and semi-analytical models is yet to 
be seen. This and similar approaches only partly address 
the issue whether $(L\M,X) = $(L|M) since when one 
redistributes galaxies to different halos one cannot study 
differences coming from the internal structure of that halo. 
For example, older halos may have a dark matter profile 
with a larger concentration parameter and galaxy proper- 
ties may reflect that compared to a same mass halo but 
formed later. When one compares a redistributed sample 
to the real one, effects coming within halos will simply be 
averaged out. Addressing whether $(L|M, X) = <&(L\M) 
is challenging, but large catalogs of galaxies from surveys 
such as SDSS and 2dF provide the opportunity for a va- 
riety of studies. Though how these catalogs will provide 
the ultimate answer is yet to be seen. 

One may be able to establish few percent differences 
with dependences in the environment through a combina- 
tion of luminosity functions and clustering as a function 
of overdensity. These statistics depend on different com- 
binations of dn(S ga \)/dM and $(L|M, <5 ga i) (see, Cooray 
2005c for expressions), for example, and the combination 
of one-point (luminosity function) and two-point (cluster- 
ing statistics) allows degeneracies to be broken. Still, the 
test is challenging. Based on a simple Fisher matrix test, 
using the concentration dependent bias and occupation 
properties of Wechslcr et al. (2005) as the basis for vari- 
ations with age, we find that one is required to have 1% 
accurate measurements in clustering statistics and lumi- 
nosity functions divided to ~ 10 bins of some parameter 
that is taken to be a direct estimator of age. The accu- 



racy of establishing differences is limited because one must 
introduce variations in bh.aio{M, 5 ga \ ore) when describing 
clustering statistics, but not in the luminosity function. 

Finally, to conclude, there is a strong possibility that 
one must account for small (at most ten percent) changes 
to galaxy properties with a second parameter in addition 
to the halo mass. This looks challenging when compared 
to the simple halo model that has a built-in assumption 
that the halo occupation distribution is defined only by 
the halo mass. If the second parameter that vary prop- 
erties of the conditional functions can be established, this 
information can be easily included when interpreting data 
under an assumed cosmological model or when attempt- 
ing to extract cosmology with galaxy clustering measure- 
ments. The conditional function approach, say compared 
to the halo occupation distribution, allows this to be eas- 
ily achieved since one can easily condition in terms of the 
mass and the second parameter as well as accounting for 
any joint distribution between mass and that parameter. 
We think that the conclusion that recent results invalidate 
analytical attempts to understand the galaxy distribution 
and use it for cosmological parameter estimates is prema- 
ture. 
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